Stochastic vortex dynamics in two-dimensional easy-plane ferromagnets: 

Multiplicative versus additive noise 
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We study how thermal fluctuations affect the dynamics of vortices in the two-dimensional clas- 
sical, ferromagnetic, anisotropic Heisenberg model depending on their additive or multiplicative 
character. Using a collective coordinate theory, we analytically show that multiplicative noise, aris- 
ing from fluctuations in the local fleld term of the Landau-Lifshitz equations, and Langevin-like 
additive noise both have the same effect on vortex dynamics (within a very plausible assumption 
consistent with the collective coordinate approach). This is a non-trivial result, as multiplicative 
and additive noises usually modify the dynamics quite differently. We also carry out numerical 
simulations of both versions of the model finding that they indeed give rise to very similar vortex 
dynamics. 
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I. INTRODUCTION 



In a large class of nonlinear problems arising in physics, 
chemistry and biology, coherent, localized excitations of- 
ten play a crucial role in governing the dynamics of 
the corresponding systems. This is the case, for in- 
stance, with solitons, vortices, fronts and many other 
solitary wave-like objects, found in a wide variety of 
low dimensional systemsEd Physical situations featur- 
ing these phenomena are usually described by one of a 
few "canonical" partial differential equations, either in- 
tegrable, such as the one-dimensional sine-Gordon or the 
nonlinear Schrodinger equations, or non-integrable, such 
as the (/)^ one.cl However, those mathematical formula- 
tions correspond in general to highly idealized physical 
contexts, in which factors like inhomogeneities, fluctu- 
ations, external fields, or damping are not taken into 
account. Conclusions about stability, dynamics, inter- 
actions and any other property of coherent excitations 
drawn from those simple descriptions do not necessarily 
carry over when the above ingredients can not be ne- 
glected. As a consequence, interest on the effect of per- 
turbations on solitons and related excitations has grown 
rapidly since the early eighties, motivated by the need of 
bridging the gap between ideal models and real applica- 
tions. 

One of the most important and universal perturba- 
tions of any physical system is noise in one of its various 
forms.Ll Typically, noise enters the physics of a system ei- 
ther as temporal fluctuations of internal degrees of free- 
dom, caused by temperature, for instance, or as random 



variations of some external control parameter. In order 
to study the effect of these fluctuations in the system one 
is interested in, random terms are added to the original 
deterministic equations; generally speaking, internal ran- 
domness will reflect itself in additive noise terms, while 
external fluctuations will almost always give rise to mul- 
tiplicative noise terms. The difference between the two 
cases is that additive noise describes fluctuations inde- 
pendent of the values of the system variables, whereas 
multiplicative noise relates to fluctuations whose magni- 
tude is modifled by the state of the system. Of course, 
adding one kind of noise or the other to an otherwise 
deterministic problem leads in general to very different 
results, and nonlinear coherent excitations are not an ex-, 
ception: Thus, for instance, studies of the sine-GordonH 
and (j)'^ (Refs. ItIJI) equations have shown that large scale 
(i.e., spatially homogeneous) noise modifles the dynamics 
of solitons in very different ways depending on its addi- 
tive or multiplicative character. As another example, we 
note that the difference between additive and multiplica- 
tive noise in the nonlinear Schrodinger equation has also 
been discussed in Ref. |9| where multiplicative noise is as- 
sociated with scattering of solitons by phonons with no 
creation of new phonons, whereas additive noise implies 
creation and anihilation of phonons. However, a gen- 
eral discussion of the analogies and differences of both 
kinds of fluctuations is lacking in nonlinear partial dif- 
ferential equations. Moreover, studies of noise effects on 
model systems are often carried out without any reference 
to the physical meaning of the type of term introduced; 
hence, we believe that a physically clear-cut example will 
help understand the similarities and differences of addi- 



tive and multiplicative noise in other cases. 

In this paper, we aim to gain insight into the effects of 
the two types of fluctuations by presenting a comparison 
of the effects of additive and multiplicative noise on the 
dynamics of vortices in two-dimensional (2D) easy-plane 
ferromagnets, as described by the classical, anisotropic 
Heisenberg model. The reason we choose this system is 
that we can justify physically in very direct ways the 
reasons for introducing one or the other type of noise in 
the Landau-Lifshitz equation, thus making possible the 
discussion and interpretation of our results in physical 
terms. To this end, we will address the problem both 
from analytical and numerical viewpoints in order to 
achieve a more complete picture of the two cases. Accord- 
ingly, in Section 2 we introduce our model, summarize its 
main features, and discuss how noise can be introduced 
in either form according to the physics one has in mind. 
In Sec. 3, we present our analytical results, obtained in 
the framework of a collective coordinate approach. This 
approach will allow us to show the surprising result that, 
with a reasonable assumption, very natural in the context 
of a collective coordinate theory, the two kinds of fluctu- 
ations considered lead to the same dynamics for a single 
vortex. Our analytical predictions are tested in Sec. 4 by 
means of numerical simulations, which yield a very good 
agreement with the theory thus confirming a posteriori 
the validity of our assumption. Finally, Sec. 5 discusses 
our results and summarizes our main conclusions. It is 
important to note that, in principle, a formulation alter- 
native to the one presented here is possible in terms of 
the Hamilton equations, instead of the Landau-Lifshitz 
equation. We discuss this possibility in an appendix and 
show that it suffers from several problems. 



II. MODEL AND STOCHASTIC 
PERTURBATIONS 

The model we will be working with is the 2D 
anisotropic Heisenberg model with XY- or easy-plane 
symmetry, given by 
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where the subindices x,y or z stand for the spin compo- 
nents, < ^ < 1, and < m,n > labels nearest neighbors 
of a square lattice. Among its excitations, specially inter- 
esting ones are vortices, that are planar (i. e., with null 
z component) if (5 > 0.297 and nopaplanar (i. e., with lo- 
calized Sz structure) ii S < 0.297OE3. Such non-planar 
vortices will be the specific object of our study as re- 
ported in the remainder of the paper; however, the ideas 
we will be discussing are general enough to be of interest 
in other, related contexts where the system behavior is 
governed by soliton-like collective excitations. 

Physically, this model has many interesting applica- 
tions: In the last few years several classes of materials 



have been found or fabricated for which magnetic inter- 
actions within planes of their crystalline structure are 
much stronger than between these planes, and therefore 
the magnetic properties are basically 2D. Materials in 
these classes include, for instance, layered magnets (such 
as Rb2CrCl4), graphite intercalated compounds (such 
as C0CI2), magnetic lipid layers (such as manganese 
stearate), and high Tc superconductors (see references in, 
e. g., Ref. |l^). It is evident that in order to model better 
these materials, one of the factors that has to be added 
to a description in terms of the Heisenberg Hamiltonian 
is fluctuations, which can arise from different origins. 
Among those, thermal noise is of course the most natu- 
ral source of fluctuations to consider in the context of the 
Heisenberg model: Indeed, from the experimental point 
of view, insofar as the motion of vortices has rueasur- 
able consequences in inelastic neutron spitteringlj and 
nuclear magnetic resonance experiments,lij the effects of 
finite temperature on vortex dynamics can have signa- 
tures in those measurements. 

In most cases, thermal effects are studied by adding to 
the equations of motion of the system under considera- 
tion an additive noise term plus a damping term in order 
to ensure that the fluctuation-dissipation theorem holds. 
This is simply the familiar Langevin dynamics. In our 
case, our starting (deterministic) point is the Landau- 
Lifshitz equation, which reads 
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where Sm is the spin vector at lattice site m, and H 
is the Hamiltonian, in our case that of the anisotropic 
Heisenberg model, Eq. (|l|). The corresponding Langevin 
dynamics equation for our model is obtained by adding 
damping and additive noise, which yields 
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The second term on the rhs of Eq. ( pi) is the damp- 
ing term. Following Refs. |5[ p^r^d p7| we have cho- 
sen for simplicity Gilbert damping ,113 chiefly because it is. 
isotropic in contrast to the Landau-Lifshitz damping.ta 
The last term in the rhs of Eq. (j^) is the noise term, given 
by a Gaussian white noise with statistics defined by 
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(4b) 



where D = 2ekBT is the diffusion constant and a, (3 = 
1, 2, 3 denote cartesian coordinates. It is important to 
note that, strictly speaking, the three Eqs. (||) do not rep- 
resent Langevin equations, because all the components of 
dSm/dt appear in each equation due to the cross prod- 
ucts. To properly introduce the noise, one first has to 
group all the time derivatives on the Ihs of the equa- 
tion, and only then one can add independent white noise 



terms, say pm, for each spin component. In Refs. ^ and 
pi| it was shown that such a procedure leads to a different 
Langevin equation, in which the noise term ?/„ must be 
replaced by fjm = Pm — ^{S x pm)- However, as discussed 
in Refs. ^ and ^, the correction is of the order of e^, 
and taking into account that in the simulations values of 
e of the order of 10"'^ are used, the correction factor can 
be neglected. I. e. we take Eq. (g) as Langevin equation, 
containing purely additive noise. 

What we have discussed is the usual way to introduce 
thermal fluctuations in any model. Nevertheless, in our 
case this Langevin approach suffers from the drawback 
that it enlarges the length of the spins, which in the orig- 
inal Heisenberg model is fixed to \S\ — 1. In the simula- 
tions of Ref. ^, this unphysical effect was suppressed by 
renormalizing the length of every spin back to its origi- 
nal length at every time step. A more detailed discussion 
of the implementation of this constraint can be found in 
Ref. g^. At this point, we were not fully satisfied with this 
solution, and looked for another way to study thermal ef- 
fects which preserve the spin length exactly. With this 
motivation, we now propose to introduce a noise term 
multiplicatively, according to 
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The term hm{t) is again a set of independent Gaussian 
white noises, but now they represent fiuctuations in the 
local field, dH/dSm, in which the spin Sm precesses. 
In fact, this is a natural way to introduce the effect of 
thermal fluctuations, as the local fleld is the only way 
through which the spin Sm can feel any changes in its en- 
vironment, due to those thermal fluctuations or for any 
other reason. Thus, the random term accounts for the 
interaction of the spin degrees of freedom with phonons, 
magnons, and any other excitation thermally generated. 
In addition, this fashion of introducing the noise has the 
property that Eq. (ph exactly preserves the spin length, 
hence there is no need for corrections as in the addi- 
tive case. We note that a similar term has been consid- 
ered by Garanin,tJ who proposed it in order to derive an 
all-temperature theory from the corresponding Fokker- 
Planck equation, obtaining a so-called Landau-Lifshitz- 
Bloch equation. Our purpose now is to use multiplicative 
noise in order to understand the influence of flnite tem- 
perature on vortex dynamics, comparing the results with 
those arising from the usual Langevin approach discussed 
above. 

There is an important question that deserves discus- 
sion before proceeding to the study of vortex dynamics, 
namely the correct way to interpret the stochastic partial 
differential equations (|^) and (H). The first one contains 
only additive noise, whick iaiplies that Ito or Stratonovic 
interpretations coincide,Qa and therefore there is no 
problem in that case. As for the second one, being multi- 
plicative, we do have to specify our interpretation of the 



equation. In principle, when thinking of thermal exci- 
tations interacting with the spins, we would have to as- 
sociate with them a finite correlation time which would 
lead to a colored noise term. Taking white noise means 
taking the limit of zero correlation time, and therefore it 
is necessary to interpret Eq. (^ in the Stratonovic sense. 
Another reason for us to stay with this interpretation is 
that the spin modulus is conserved; it can be seen that 
Ito calculus leads to an exponential decrease of the mod- 
ulus, with a damping time proportional to the damping e 
in Eq. m^In addition, in the Stratoapvic intese^retation, 



Gara 



and even earlier Brown,cil showedE3 that the 



stationary solution to the Fokker-Planck equation corre- 
sponding to Eq. (ph is the Boltzmann factor with Hamil- 
tonian (|l|), indicating that Eq. (||) indeed represents the 
dynamics of our model at finite temperatures. Interest- 
ingly, it is not difficult to show that the Boltzmann factor 
is not a stationary solution to the Fokker-Planck equation 
for the additive noise Langevin problem, Eq. (^, unless 
further assumptions are |tqade, including the constraint 
of constant spin modulu^j. 



III. ANALYTICAL RESULTS 

As stated in the Introduction, our approach to the 
problem of vortex dynamics will be both analytical and 
numerical. In this section, we first derive equations of 
motion for the vortex center X{t) both for the additive 
and for the multiplicative noise cases, and afterwards we 
compare with numerical simulations for our model, i. e., 
with results from numerical integration of Eqs. (0) and 

i). 

Our analytical approach to the stochastic dynamics of 
vortices begins by taking the continuum limit of Eqs. 
(g|) and (|5|), which are much more difficult to deal with 
in a lattice formulation. This is a good approximation 
provided that the localized S'^-structure spans many lat- 
tice sites, as typically occurs in practice (except if the 
anisotropy parameter is chosen close to the critical value 
S — 0.297). The next step is to use a collective co- 
ordinate theory to analyze the vortex dynamics, (see 
Ref. p^ for a recent review on collective coordinate ap- 
proaches). Within this procedure, one assumes that 
the shape of the excitation under consideration, in our 
case a vortex, is not modifled by the perturbation for a 
large range of pectiirbation types (this is a very general, 
widely applicableEZI approach), i.e., in our case noise and 
damping terms, and that only the dynamics of its cen- 
ter is modified by these extra terms. The vortex mo- 
tion is then introduced by the travelling wave Ansatz 
S{f,t) — S{r ~ A(t)), where S{r) describes the static 
vortex shape. Unfortunately, such a simple apopaach 
(first proposed for magnetic domains by Thielell3't3) is 
not enough to describe the vortex dynamics, as was found 
in Ref. pSl. In that work, Mertens et al. developed a gen- 
eralization of the collective coordinate theory in which 



the vortex shape is allowed to depend on the velocity X 
and, in general, also on higher order derivatives of X{t). 
The corresponding generalized travelling wave Ansatz is 



5(r,t) = 5(r-X,X,X,...,X("^), 



(6) 



which yields an [n + l)-th order differential equation for 
X{t). However, as discussed in Ref. ^, in the case of 
non-planar vortices only the odd-order equations repre- 
sent self-consistent valid approximations; and it turns out 
that the third-order equation is sufficient to describe ac- 
curately all simulations without dampingc3. Therefore, 
in this paper we use the Ansatz ffl) with n — 2 and apply 
it to the general case, which includes damping and noise. 
The continuum versions of Eqs. (0) and (0) read, re- 
spectively 
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To obtain the equations for the collective coordinate, in- 
stead of using the Hamiltonian procedure described in 
Ref. 28, we follow a much more direct approach, which 
we have already used for the additive noise (see a prelim- 
inary report in Refs. ^ and ^9|): We begin with Eq. (0) 
and multiply it by S" ■ (dS/dXi) x , where Xi is the i-th 
component of the vortex center position. The contribu- 
tions of all terms in the rhs of Eq. (q) are 
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where Ti is the Hamiltonian density in the continuum 
limit of Eq. (nl). According to our ansatz we insert in Eq. 
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The Ihs of Eq. (|^) is dealt with in the same way as just 
described for the rhs. By collecting the results, inte- 
grating over r, and dividing by S*^, we obtain the same 
third-order equation as that in Refs. ^,^, and |29|: 



(A + a) X + (M + m) X + (G + g) X = 

= kX + MX + QX = F + F 



add 
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The terms in Eq. (11) are as follows: The tensors denoted 
by capital letters come from the Ihs of Eq. (0) , and their 
expressions are, for the gyrotensor G, 



Gy- = S-" I d'rS- 



for the mass tensor M, 
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and for the third-order gyrotensor A, 
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The tensors denoted by small letters come from the 
Gilbert damping term; as can be seen from Eq. (11), 
they contribute to all orders, and they are given by 



Qi-j = e i d^r 



dS dS 



dX, dX, 



= e Id^r 



dS dS 



dX^ dX,' 



, ,2 dS dS 



dX, dXj 



Finally, the force terms are 
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The key to achieve a complete understanding of the vor- 
tex dynamics as described by Eq. (11), is to know the 
mean (Ff'*'^) and the variance Var(Ff'^'^). The mean is 
easily shown to be zero, whereas for the correlation func- 
tions, by using the continuum version of Eq. (4b) we 
obtain 

We now turn to the multiphcative noise case, Eq. (||). 
The calculations are mostly the same as described above, 



except that the contribution in Eq. (^) is now substi- 
tuted by the term coming from the multiphcative noise, 
which reads 
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As a consequence, only the stochastic term in Eq. (11) is 
modified: The new stochastic force is 
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which is in principle different from what we obtained for 
the additive case, F^'^'^. However, when we evaluate the 
first moments, we find again that the mean of i^™"!* is 
zero. Furthermore, to compute the variance, we need to 
evaluate 
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where summation over a and j3 is implied. At this point, 
it is important to notice that the noise is multiplicative, 
and therefore in principle we cannot take the spin fields 
out of the average over realizations of the noise. We stress 
that this is not a problem in the case of additive noise, 
and that we indeed proceeded that way to obtain the ex- 
pression for the variance of F^'^'^ given in Eq. (20). This 
is so because the additive character of the noise implies 
that the spin fields and the noise are uncorrelated. How- 
ever, one can not simply apply the same argument to the 
calculation in Eq. (23), because when the noise is mul- 
tiplicative it is not clear whether the spin fields and the 
noise are correlated or not. In this situation, in principle 
we can not exactly evaluate the variance of F™^'*^ but 
we can make the following approximation: At least for 
small noise, we can substitute the spin fields in Eq. (23) 
by the deterministic expression for the vortices, assum- 
ing reasonably that the corrections induced by the noise 
will be of the order of the noise strength and that their 
contribution to the variance would be two orders higher 
in the noise strength than that of the deterministic part, 
and hence, negligible. If we do so, we can then take the 
spin fields out of the average, finding 
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i.e., the variance of the stochastic force due to the multi- 
plicative noise is the same as that of the additive noise. 
We stress that the assumption leading to the above, 
unexpected result, is very natural within a collective co- 
ordinate approach such as the one we are using here. 
To understand this, recall that the main hypothesis of a 



collective coordinate theory is that the shape of the exci- 
tation under study remains mostly unaffected by the per- 
turbations, and only its position and possibly a few other 
collective variables change due to therriEj This assump- 
tion amounts to a drastic reduction in the degrees of free- 
dom of the system: From an infinite number of them in 
the continuum equations to two for the center dynamics, 
in the present case. Physically, ii-is equivalent to neglect 
the contribution of the magnona^J excited by the pertur- 
bations to the vortex shape; and this approximation is 
made already at the very beginning, in order to obtain 
Eq. (11). In view of this, it would not be reasonable 
to keep the exact fields S, containing the contribution of 
the phonon degrees of freedom, in Eq. (23), and therefore 
we carry out the calculation with the deterministic vor- 
tex shape, obtaining the result (24). Ell Thus, consistently 
within our collective coordinate approach, we have shown 
that the mean and the variance of the stochastic forces 
induced by the multiplicative and the additive noises are 
the same (higher-order moments may still differ). As a 
consequence, the dynamics of a vortex is predicted to be 
the same under the influence of each type of noise, which 
is a surprising result in view of the general result that 
their effects are very different. 

Once we have shown that both stochastic forces are 
equal, then Eq. (11) has to be the description of the vor- 
tex dynamics under the two types of noise. This equation 
was already studied and solved in Refs. ^ and ^, and we 
will only summarize here what we need for our present 
purposes. The classical spin is constrained to have a fixed 
magnitude which we set to unity. Therefore, the integrals 
on the rhs of Eqs. (|2|) through (|7|) were calculated us- 
ing canonical fields = arctan(S'y/S'a;) and ip = Sz for 
the spin vector: 



S = v 1 — "0^ cos (f) Cx + \/ 1 — ijj'^ sin (f>ey -h ^JCz- (25) 

Now, the explicit calculation of all the integrals is only 
possible if the dynamic structure of the vortex is known. 
In this respect, in Ref. |2^ it was shown that the core re- 
gion of the vortex contributes very little [except to Eq. 
([l2|)]; the dominant contributions stem from the outer 
region, if the system size is large enough. A vortex in 
the center of a circular system with radius L and free 
boundary conditions has the following structurepin the 
outer region which was confirmed by simulations n3 



= 4>0 + (f>l + (j)2, -0 = V'o + V-'l + "02 (26) 



with 



A ^ -1 3^2 

bo = q tan — , 

Xi 



(fii ^p{xiXi +a;2X2), 



q r •• 

^2 = -BT log ^-(a;2Xi - 0:1X2) 
80 eL 



(27) 



(28) 



(29) 



and 



V'o 



ipi 



■02 



P\I — '^^P( 



-r/ry), 



4(5r- 



r(x2Xi - X1X2) 



4(5 



(xiXi + 2:2X2 



(30) 



(31) 



(32) 



Here q = ±1 is the vorticity and p — ±1 is the po- 
larization, which determines to which side the out-of- 
plane structure of the vortex points. In addition, r^, — 
[(1 — (5)/5]^/^/2 characterizes the radius of the vortex 
core. Straightforward integrations then yield expressions 
for the tensors G, M, and A, as well as for the forces, 
which can afterwards be inserted in Eq. (11) for the vor- 
tex motion. Finally, this equation is linear except for the 
force F(X), which can in turn be linearized by expanding 
around the mean trajectory. Subsequently, the equation 
can be solved by means of a Green's function approach 
and the so obtained solution can be used to calculate 
analytically the variances of the vortex trajectory X(i). 
These variances are proportional to the effective vortex 
diffusion constant 



Dy^Diri log — + Ciac) 
an 
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where ac is chosen of the order of r^, , and the unknown 
constant C stems from the core region. As in this work 
our main interest is to show that the dynamics under ad- 
ditive and multiplicative noises is the same, we will not 
dwell any further in the solution of Eq. (11); the inter- 
ested reader is referred to Refs. ^ and ^ for details. 



IV. NUMERICAL RESULTS 



our interpretation of all the equations in the sense of 
Stratpnoyic, we use the Heun method to integrate in 
timcaa We explicitly take into account the constraint 
5^ = 1 by means of a Lagrange parameter, as discussed 
in Ref. ^3 (see also Ref. Elf). To find a proper damping 
constant we checked the time dependence of the system 
energy using different damping constants for L — 2A and 
T — 0.02 (in dimcnsionlcss units). The energy at t = is 
the same as for T = and e = because the noise is in- 
troduced with the first time step of the simulation. The 
energy then rises and saturates to a value independent of 
e, but for e > 8 X 10^"^ the energy decreases slowly after 
saturation. The saturation time gets longer with lower e, 
for e > 2 X 10^'^ we achieve acceptable saturation times 
< 300 (in units of h/{JS)). We have always made a 
pre-run of length to > 300 prior to beginning the eval- 
uation of the simulation data. Finally, our simulations 
consisted of numerical integrations up to times f = 4 000 
(note in this regard that this takes 10 days CPU time on 
a 433-MHz-Digital-Alpha workstation for averages over 
100 runs) because this is larger than the characteristic 
time given by 5/e for the damping in the trajectories. Fi- 
nally, the difference between the energy without temper- 
ature and the saturation energy with temperature must 
be the thermal energy. We computed the mean thermal 
energy per spin at several temperatures and it agreed 
with //2 X ksT up to T = 0.9 (for comparison, we note 
that the Kosterlitz-Thouless transistion temperature is 
about 0.8 in our units), / being the number of degrees of 
freedom per spin. This is far above the temperatures we 
will discuss below. We have to mention that for T > 0.1 
there is a more complex phenomenology because the ther- 
mal noise does not only induce a diffusive motion of the 
vortices but it can also ffip their out-of-plane structure, 
correspondingly changing their direction of motion, and 
also nucleate additional vortex-antivortex pairs. In this 
work we stay away from this regime which has been con- 
sidered in Ref. |3^. 



Our numerical simulations begin with one vortex with 
its center located at a distance Rq from the middle of a 
circularly shaped square lattice with a radius of L lattice 
constants. We use free boundary conditions to produce 
an image ^ajitivortex which leads to a radial force on our 



vortexEj'Lj The initial spin configuration is determined 
from an iterative program which produces a discrete vor- 
tex structure on the lattice.tS In this way we avoid the 
radiation of spin waves which would appear during the 
early time units if we use a continuum or other approx- 
imation for the vortex shape. The parameters we used 
in our simulations were: S = 0.1, for the radius of the 
out-of-plane structure of the vortex to be large enough 
to avoid severe discreteness effects; system radius L = 24 
for the vortex to have enough space to evolve far from 
the boundaries; and initial distance from the circle center 
Ro = 10. 

For the time integration of the Landau-Lifshitz equa- 
tion we use the discrete version of (y) . Consistently with 



The outcome of our numerical simulations can be sum- 
marized by saying that they fully confirm the predictions 
of our analytical calculations, namely that the vortex 
dynamics is the same under both kinds of noise. An 
example of our results is shown in Fig. yj for additive 
noise and in Fig. for multiplicative noise. It is already 
evident from comparing the dispersion in the R and (j) 
components that the behavior of the vortex is very ap- 
proximately equal in the two problems. 




FIG. 1. Variances of the vortex trajectory for Langevin dy- 
namics with additive noise. The temperature is T = 0.03; 
other parameters are given in the text. From top to 
bottom, shown are the variance of the radial coordinate, 
°"fli? = (-R'^) ~ {R^y '^fl0i the off-diagonal elements of the 
variance matrix, and o-^j,, the variance of the azimuthal co- 
ordinate. In all three cases the lower line is the theoretical 
prediction with the vortex diffusion constant chosen in order 
to make the analytical curve lie just below the simulated data. 




FIG. 2. Variances of the vortex trajectory for the case of 
multiplicative noise; parameters are the same as in Fig. y. 
From top to bottom, shown are the variance of the radial co- 
ordinate, (jjiji — (i?^) — (R)'^, '^R</>i the off-diagonal elements 
of the variance matrix, and a^^, the variance of the azimuthal 
coordinate. In all three cases the lower line is the theoreti- 
cal prediction with the vortex diffusion constant fitted to the 
data. 



Furthermore, the comparison to the analytical predic- 
tion is very good: We emphasize that there is only one ad- 
justable parameter, the vortex diffusion constant, whose 
value we can only estimate as an exact expression for the 
discrete structure of the vortex core is not known. The 
systematic deviation of the theory_for very long times, 
t > 2000, is due to a simplification,Ej namely keeping Rq 
constant; in the simulation, Rq slowly increases because 
of the damping. Results for other low temperatures com- 
pare equally well, the better the lower the temperature; 
for higher temperatures, our numerical estimates are less 
accurate, although qualitatively the results for both cases 
remain the same. 



V. CONCLUSIONS 

In this work, we have provided analytical and numeri- 
cal evidence that additive (Langevin-like) and multiplica- 
tive noise (coming from fluctuations in the local field) 
have the same effects on the dynamics of vortices de- 
scribed by the 2D Landau-Lifshitz-Gilbert equation. An- 
alytically, the result has been obtained in the framework 
of a collective coordinajte. approach within a generalized 
travelling wave AnsatzEB The variances of the effective 
force acting on the vortex were shown to be the same in 
the additive and multiplicative cases provided that de- 
formations of the unperturbed vortex shape can be ne- 



glected. It is important to stress that this hypothesis is 
not an extra assumption but rather it is in fact imphcitij^ 
made when using any collective coordinate approach.EZI 
In order to substantiate our analytical results, we numer- 
ically simulated both the additive and the multiplicative 
cases for the 2D easy-plane Heisenberg ferromagnet, find- 
ing an excellent agreement with the prediction of equal 
behavior under both sources of noise, as well as with the 
analytical expression for the variance of the vortex tra- 
jectory. 

As an immediate consequence of the validity of the 
above discussed prediction, we point out that all the re- 
sults obtained in Ref. ^ for the additive (thermal) noise 
apply to the multiplicative model presented here, specif- 
ically: i) the existence of three different temperature 
regimes for the vortex propagation: a low temperature 
one, where the vortex motion follows essentially the third 
order equation of motion with parameters independent of 
temperature; a middle temperature one, at which traces 
of the oscillations arising from the third order equation 
are lost, and a high temperature regime, which is not 
describable by a one-vortex approach because too many 
vortex-antivortex pairs arise in the system; and ii) the 
dependence of the effective diffusion coefficient for the 
vortex on temperature. On the other hand, the work we 
have reported on here allows to place on firmer grounds 
that those are indeed the features of thermal vortex dy- 
namics: The problem with the non-constant spin length 
in a Langevin approach is now solved by the multiplica- 
tive approach, in which it is exactly conserved. Further- 
more, the fact that the stationary probability distribu- 
tion for the multiplicative case is precisely the Boltzmann 
factor reinforces our conclusion that this is the correct 
description of thermal effects in the framework of mod- 
els with dynamics given by the Landau-Lifshitz-Gilbert 
equation. The thorough study carried out in Ref. ^ for 
the Langevin approach is then fully confirmed by the 
present work. 



ACKNOWLEDGMENTS 

We thank Grant Lythe, Yuri Gaididei, Francisco 
Domingucz-Adame, Raiil Toral, and Niels Gr0nbech- 
Jenscn for discussions. 

Travel between Bayreuth and Madrid is supported by 
"Acciones Integradas Hispano-Alemanas" , a joint pro- 
gram of DAAD (Az. 314-AI) and DGES. Additional 
funding for the stay of E. M. at Bayreuth was provided 
by the Fundacion Universidad Carlos III. Travel between 
Europe and Los Alamos is supported by NATO grant 
CRG 971090. Work at Leganes is supported by CICyT 
(Spain) grant MAT95-0325 and by DGES (Spain) grant 
PB96-0119. Work at Los Alamos is supported by the 
United States Department of Energy. 



APPENDIX A: NOISE IN THE HAMILTON 
EQUATIONS 

In this section we consider an alternative formulation 
of the problem of thermal fluctuations in the classical 
Heisenberg ferromagnet, and show that this version ex- 
hibits different features for additive and multiplicative 
noises, not even being well posed for the latter case. 

By taking into account that the spin length has to be 
constant, one can reformulate the Heisenberg model in 
terms of the fields ((> = arctan(S'j,/S'a;) and ijj = Sz for the 
spin vector: In this way, after rewriting the Hamiltonian 
(|l|) in these new variables, the dynamics of the model is 
governed by the following Hamilton equations: 



t/j 



5H_ 
6H 



(Ala) 
(Alb) 



The additive noise version of the problem, i.e., the 
Langevin equations, is given by 



5H 
5H 



V<pir,t), 



■0 = ^"^+'?'A(^'^)' 



whereas the multiplicative version corresponds to 
SH 



Tp = -^ +'l]^{f,t)'lp. 



(A2a) 
(A2b) 

(A3a) 
(A3b) 



Of course, both models should also be supplemented 
with damping terms in order to fulfill the fluctuation- 
dissipation theorem. For simplicity, we will leave them 
out in this discussion. Their influence, however, is not 
qualitatively important as for our main argument below, 
as it can be checked that they would only contribute fac- 
tors to the expressions for the force variances. 

Now, from either one of these two equations we can 
derive an equation of motion for the vortex center fol- 
lowing basically the same collective coordinate approach 
as above. In this case, one has to multiply the equation 
for ip by d4>/dXi and substract from it the equation for 
(j) times dip/dXi. In this way an equation completely 
analogous to the third order equation of motion (11) is 
obtained, but now we flnd for the forces in the additive 
case 



F 



add 



- d'r 



dX, 



■m 



dip 
dxi 



m 



(A4) 



whereas for the multiplicative noise. 



pnmlt ^ / ^2^ 



d(j} ^V' , ' 



(A5) 



Let us first discuss the additive case, Eq. (A4). It 
is not difficult to derive from the above expression the 
first moments of the additive force, which turn out to be 
{F^'^'^it)) =Oand 



DS{t - 1') jd? 



dX,, 



f dip 
\dX, 



(A6) 



This expression is actually not very different from that 
obtained in Eq. (20) for the additive noise version of the 
Landau-Lifshitz equation: Using Eq. (P5|), the integrand 
of Eq. (20) reads, for i — j, 



m ^<-^Mt 



1 



dip 



(A7) 



As the static Sz structure of the vortex (^0|) falls off ex- 
ponentially in the outer region Uc < r < L, where L is 
the system radius, both integrals (20) and (A6) have the 
same logarithmic size dependence, 

^^add^^^^add^^/^^ = D5{t - t')7T flog —+C 



(A8) 

The only difference between both formulae consists in 
the value of the (small) constant C, which stems from 
the core region < r < Oc, where Oc = 0{ry). 

M oving now to the multiplicative case, given by Eq. 
(A5), it can be easily shown that the situation is com- 



pletely different. The mean and the variance of F™" can 
again be easily computed, yielding (i^™""(0) — and 

(i^f""(i)i^f""(t')) = 
D6it - 1') I 



(fr 



\dX, 



^' 



dip 



(A9) 



which, due to the presence of a factor ip'^ in the first term 
of the integrand decaying exponentially away from the 
vortex core, is independent of the system size. It is then 
clear that the variances are very different in both cases. 
Therefore, a detailed stud y of the Hamilton equations 
with additive noise, Eqs. (A2) would basically lead to 
the same conclusions summarized above arising from the 
study of the Landau-Lifshitz-Gilbert problem, whereas 
the multiplicative version would give totally different an- 
swers. It has to be concluded then that the proper way 
to include multiplicative noise in this problem is at the 
Landau-Lifshitz level and not in the Hamilton equations 
for (j) and ip. Finally, this result points out very clearly 
that the fact that the two approaches, additive and mul- 
tiplicative, give the same results in the Landau-Lifshitz 
equation, is neither trivial nor generally applicable, and 
should then be regarded as a specific and attractive fea- 
ture of that approach. 
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